Helical filaments with varying cross section radius 
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Abstract 

The tridimensional configuration and the twist density of helical rods with varying cross section 
radius are studied within the framework of the Kirchhoff rod model. It is shown that the twist 
density increases when the cross section radius decreases. Some tridimensional configurations of 
helix-like rods are displayed showing the effects of the nonhomogeneity considered here. Since the 
helix-like solutions of the nonhomogeneous rods do not present constant curvature and torsion a 
set of differential equations for these quantities is presented. We discuss the results and possible 
consequences. 
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Helical filaments are tridimensional structures universally found in Nature. They can 

n q 

be seen in very small sized systems, as biomolecules [ll and bacterial fibers |2j, and in 
macroscopic ones, as ropes, strings, and climbing plants J, U, All these objects have in 
common the fact that the mathematical geometric properties of the 3D-space curve related 
to their axis, namely the curvature, kp, and the torsion, Tp, are constant jf|. 

The so-called Kirchhoff rod model has been proved to be a good framework to study the 
statics laLDl and dynamics M of long, thin and inextensible elastic rods 9UlC|. The applica- 
tions of the Kirchhoff model range from Biology pj, |5|, [11|] to Engineering 12J. In these cases, 
the rod or filament is considered as being homogeneous, but the case of nonhomogeneous 
rods have been also considered in the literature. It has been shown that nonhomogeneous 
Kirchhon rod8 may present spatial chaos Q Q and that helical transitions occn r in the 
tridimensional configurations of rods with periodic variation of the Young's modulus Q|. 
A comparison between homogeneous and nonhomogeneous rods subject to given boun dary 
conditions and mechanical parameters was performed by da Fonseca and de Aguiar in 
The effects of a nonhomogeneous mass distribution in the dynamics of unstable closed rods 
have been analyzed by Fonseca and de Aguiar Goriely and McMillen ^1 studied the 
dynamics of cracking whips and Kashimoto and Shiraishi studied twisting waves in 
inhomogeneous rods. 

Here, using the Kirchhoff model, we shall present the results for the equilibrium solutions 
of nonhomogeneous rods with varying cross section radius and no intrinsic curvature. Only 
the solutions classified by Nizette and Goriely [f| as being helical will be considered: the 
straight rod, the twisted planar ring and the helix. We shall show that the twist density 
varies along the rod inversely proportional to the fourth power of the radius of the cross 
section. Also, it will be seen that the curvature, kp, and the torsion, t f , are not constant 
for the helix-like solutions with nonhomogeneous cross section radius. 
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The motivations for this work are: i) the study of failure or rupture of cables 
For example, it was shown that shoreline anchor rods rupture at the region where the rod 
diameter diminishes due to corrosion j^]. The fact that, for twisted rods, the twist density 
increases in the regions where the rod diameter decreases can be related to the onset of the 
failure; ii) the shape of some climbing plants have filamentary helical structures (spring-like 
tendrils) whose radius and pitch are not constant ^3]. Such a tridimensional configuration, 
with the radius and the pitch varying along the rod, will be shown to be a possible solution 



of the Kirchhoff model for the nonhomogeneous rod. Other motivations are related to 



defects |23j], distortions [24( and the rule of twisting 2^ in biological molecules. 

The static Kirchhoff equations, in scaled variables, for rods with circular cross section 
and no intrinsic curvature are given by: 

F' = , 

M' + d 3 x F = , (1) 
M = I(s) h di + JO) k 2 d 2 + T I(s) k 3 d 3 , 

where s is the arc-length of the rod, the prime ' denotes differentiation with respect to s 
and the vectors F and M are the resultant force, and corresponding moment with respect to 
the axis of the rod, respectively, at each cross section, dj, i = 1,2,3, compose the director 
basis with d 3 chosen to be the vector tangent to the axis of the rod and di and d 2 lie in 
the plane of the cross section, ki are the components of the twist vector, k, that controls 
the variations of the director basis along the rod through the relation d- = k x dj. k\ and 
k 2 are related to the curvature [k F = \fk\ + fcf) and k 3 is the twist density of the rod. 
T = 2fi/E is the adimensional elastic parameter, with [i and E being the shear and the 
Young's moduli, respectively. I(s) is the variable moment of inertia that is related to the 
radius of the cross section through the relation I(s) = R 4 {s) (valid in scaled units). Writing 
the resultant force F in the director basis, F = /xdx + f 2 d 2 + / 3 d 3 , the equations (JH) give 
six differential equations for the components of the resultant force and twist vector: 

f[-f2h + f 3 k 2 = 0, (2a) 

f 2 + fih-f 3 h = 0, (2b) 

fi-fik 2 + f 2 k 1 = 0, (2c) 

(I(s) h)' + (r - 1) I(s) k 2 k 3 -f 2 = 0, (2d) 

(I(s) k 2 )' - (r - 1) I(s) h k 3 + A = , (2e) 

(r i( s ) k 3 y = o . (2f) 

Since I(s) = R 4 (s), the equation (J2f|) shows that the twist density k 3 is inversely proportional 
to R A (s). Therefore, k 3 is not constant for nonhomogeneous cases. On the other hand, the 
component M 3 = T I k 3 of the moment in the director basis (also called torsional moment), 
is a constant along the rod. 



In order to look for helical solutions of the eqs. (J2J) the components of the twist vector k 
are expressed as follows: 

k\ = kp sin £ , k 2 = kp cos £ , k 3 = + t f , (3) 

where kp and Tp are the curvature and torsion of the rod, respectively. In the homogeneous 
case kp and r F are constant, and £ = (k 3 — Tp)s 

We shall consider the following cases: the straight rod, the twisted planar ring and the 
general helix-like rod. 

i) The straight rod: kp = t f = 0. 

From eq. ®, k± = k 2 = and from eqs. (J2J), 

u / \ _ Constant _ m 3 

3l6J /(,) rjt*(») , ^ 

fi — /2 = and f 3 = T — Constant . 

T is the tension applied to the rod. Figure shows the twist density k 3 (s) for a straight rod 
with the cross section radius varying as 

R(s) = 1 + 0.1 cos(0.3s) , (5) 

in scaled units. The mechanical parameters used to obtain the rod displayed in Figure Q 
were M 3 = 1CT 1 (scaled units) and T = 0.9. 

ii) The twisted planar ring: kp = Constant and Tp = 0. 

The components of the twist vector for the twisted planar ring are: 

k x = k F sm£ , k 2 = k F cos£ , f = k 3 = j^ff^y . (6) 

Substituting the eqs. (0) in eqs. (j2J shows that the twisted planar ring is a possible 
equilibrium solution only if the cross section radius is of the form: 

R{s) = (Aq cos{k F s) + B sm(k F s) + C^kp)^ , (7) 

where A , B and Cj are constants. 

Remark: considering t f = and assuming that kp is a function of s (instead of being a 
constant) there exist no solutions for eqs. (J2J). So, the existence of planar solution requires 
k F = Constant. 

hi) The general helix-like rod: kp = k F (s) and r F = t f (s). 



In this case, the eqs. (jSJ) become: 

ki = k F (s) sin£ , k 2 = k F (s) cos£ , k 3 = £' + r F (s) . (8) 

Substituting eq. (JBJ) in eqs. (J2J), extracting fi and / 2 from eqs. Pe|) and (j23J), respectively, 
differentiating them with respect to s and substituting in eqs. (|2*a|) and (J2b|) gives a set of 
differential equations for k F (s), r F (s) and / 3 (s): 

[(J(s) k F {s)){T k 3 (s) - 7>(s))]' - (I{s) k F (s))' r F (s) = , 

(I(s) k F (s))" + I(s) k F (s) t f (s)(F k 3 (s) - r F (s)) - f 3 {s) k F (s) = , (9) 
(I(s)k F (s)yk F (s)+Ms) = 0. 

These differential equations are nonlinear and depend on R(s) through I(s). 

Figure El shows a helix-like rod for the following linear variation of the radius of the cross 
section: 

R(s) = 1 + 0.0023 s . (10) 

The mechanical parameters, in scaled units, are M3 = 0.05, T = 0.9, k F (0) = 0.05, t f (0) = 
0.24 and / 3 (0) = 0. 

In the figure El we display /^(s) (full line) and 0.1 R(s) (dashed line). We can see that the 
radius and pitch of the helix-like tridimensional configuration displayed on the left of figure 
121 are not constant. 

Figure 01 shows the numerical solution for the curvature k F (s) (full line) and torsion 
t f (s) (dashed line) for the helix- like rod shown in the figure 121 Despite the simplicity of its 
tridimensional shape (figure on the left) k F (s) and r F (s) are not simple functions of the 
arclength s, showing the nonlinear characteristic of the system. 

Figure |U shows an example of a helix-like rod with periodic variation of the radius of the 
cross section: 

R(s) = 1 + 0.1 sin(O.ls) . (11) 

The mechanical parameters, in scaled units, are M 3 = 0.04, T = 0.9, k F (0) = 0.19, t f (0) = 
0.04 and / 3 (0) = 0.005. 

The functions /c 3 (s) and 0.1 R(s) are displayed in the figure |U (full line and dashed line, 
respectively). In this case of periodic variation of the radius of the cross section, the tridi- 
mensional helix-like rod displayed in figure 0] (left) is more complicated than that displayed 
in the figure 121 (left). 



Figure El shows the numerical solution for the curvature k F (s) (full line) and torsion tf{s) 
(dashed line) for the helix- like rod shown in the figure 0] The curvature and torsion are not 
simple functions of the arclength s of the rod. 

There is a kind of solution called free standing helix that is defined by setting the resultant 
force F = 2f|. In eqs. (JOJ) fs(s) = gives the following solution for the curvature and 
torsion of the rod: 



(I(s)k F (s)Y = =► k F (s) = ^, 

T F ( S )=Tk 3 (s) => T F (8) = jft. 

Notice that for the free standing helix-like rod the curvature k F (s), and the torsion t f (s), 
will be analytical functions of the arclength s if the moment of inertia J(s) is given by an 
analytic function of s. Also is a constant for all s. 

The variation of the twist density along the rod can be a key factor in a variety of 
phenomena. As mentioned before, the onset of a failure in a twisted cable can be related 
to the increasing of the twist density in a given region of the cable. In the Kirchhoff model, 
the torsional moment M 3 is constant along the rod. For a relatively high value of M 3 the 
twist density at a region of the rod with small diameter can be so large that it may not 
be valid the assumption of linear relationship between the torque and the components of 
the twist vector. Also, depending on how large the moment is, the behavior of the material 
could not be approximately elastic in the regions of small diameter. So, the increase of the 
twist density due to the decreasing diameter in a twisted rod may be the starting point of 
a process that can culminate with its rupture. 

An interesting question arises for the important phenomenon known as writhing instabil- 
ity. In this phenomenon a local change in the twist can lead to a global reconfiguration of 
the rod that is a consequence of a topological constraint given by a mathematical theorem 
by White j^. If the twist density varies along the rod, the question is to identify the region 
of the rod where this kind of instability will occur. The nonhomogeneity considered here 
can also have important consequences in the dynamics of this phenomenon [28j . It will be 
a subject of a future publication. 

Another implication of variable twist density along twisted rods is the problem of sta- 
bility of equilibrium solutions. It is known that above a critical value of the twist density 
an equilibrium solution for the Kirchhoff model becomes unstable j^, Another good 
question is to investigate if a local increasing of the twist density above the critical value, 



can lead to a global instability of the related equilibrium solution. It could be important to 
the problem of failure mentioned before. 

A very hard problem in differential geometry is obtaining a direct relationship between 
the curvature and the torsion with the radius and the pitch of a helix-like nonhomogeneous 
rod. Since the definitions of curvature and torsion involve the calculation of the modulus 
of the derivatives of tangent and normal vectors with respect to the arclength of the rod, 
for non constant radius and pitch, this relation is very complicated. The analysis of this 
problem will be considered in a future work. 

It is interesting to note that the tridimensional configuration of the figure El displays 
a pattern in the radius and the pitch of the helix-like rod seen in the spring-like tendrils 
of some climbing plants. Since the young parts of the filament that composes the plant 
has smaller diameter than the older parts, these filaments are examples of rods with the 
nonhomogeneity considered here. 

The numerical solutions obtained for Tf(s) and the solutions for ks(s) of the helix-like 
cases show that the term of the equation (jHJ) is not null as it was proved to be for the 
case of homogeneous helix It means that helix-like filaments formed by nonhomogeneous 
rods are not twistless. 

The existence of intrinsic curvature may lead to other planar solutions. The helix-like 
solutions can also be influenced by the intrinsic curvature. This will be considered in a more 
complete work. 
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FIG. 1: Top: a twisted straight rod with varying radius of the cross section. Bottom: the twist 
density as function of the arc- length s. 
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FIG. 2: Left: helix-like rod having the radius of the cross section varying linearly with s (eq. (|lUjO . 
Right: 0.1i?(s) (dashed line) and ^3(5) (full line). 




FIG. 3: Geometric parameters: curvature, kp(s) (full line), and torsion, Tp(s) (dashed line), for 
the equilibrium solution shown in the figure |2J 




FIG. 4: Left: helix-like rod having the radius of the cross section varying periodically with s (eq. 
PI). Right: 0.1 R(s) (dashed line) and k 3 (s) (full line). 




FIG. 5: Geometric parameters: curvature, /cp(s) (full line), and torsion, Tp(s) (dashed line), for 
the equilibrium solution shown in the figure |1J 



